Coherent Backscattering with Nonhnear Atomic Scatterers 



O 

o 

> 
O 

(N 
(N 



T. Wellens,^'^ B. Gremaud,^ D. Delande,^ and C. Miniatura^ 

^ Institut Non Lineaire de Nice, UMR 6618, 1361 route des Lucioles, F-06560 Valbonne 
^ Laboratoire Kastler Brossel, Universite Pierre et Marie Curie, 4 Place Jussieu, F-75005 Paris 

(Dated: February 2, 2008) 

We study coherent backscattering of a quasi-monochromatic laser by a dilute gas of cold two- 
level atoms. We consider the perturbative regime of weak intensities, where nonlinear effects arise 
from inelastic two-photon scattering processes. Here, coherent backscattering can be formed by 
interference between three different scattering amplitudes. Consequently, if elastically scattered 
photons are filtered out from the photodetection signal by means of suitable frequency-selective 
detection, we find the nonlinear backscattering enhancement factor to exceed the linear barrier two. 

PACS numbers: 42.25.Dd, 42.65-k, 32.80-t 
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Light transport inside a nonlinear medium gives rise 
to a wide variety of phenomena, like pattern forma- 
tion, four waves mixing, self focusing effects, dynami- 
cal instabilities... P]. These effects are well described and 
understood with the help of an intensity dependent sus- 
ceptibility (e.g. x^^"* nonlinearity) . However, in these 
approaches, one usually discards the fact that interfer- 
ence phenomena in disordered systems may significantly 
alter wave transport properties. From this point of view, 
the most striking systems for which one must combine 
both nonlinear and disordered descriptions are coherent 
random lasers ■ Even if in this case one would require 
an active {i.e. amplifying) medium, the key point is the 
understanding of the mutual effects between multiple in- 
terferences and nonlinear scattering. In particular, in 
the case of a dilute medium, beyond deviations from the 
linear regime for the propagation parameters (mean free 
paths, refraction indices), one may expect modifications 
of interferential corrections to transport (weak and strong 
localization regimes). A paradigmatic example is given 
by coherent backscattering (CBS) experiments, where an 
enhancement of the average intensity scattered around 
the direction opposite to the incident wave is observed 
In the linear scattering regime (where the intensity 
of the scattered wave is proportional to the incident in- 
tensity), CBS arises from interference between two re- 
versed scattering paths visiting the same scatterers, but 
in opposite order. In this case, the CBS enhancement 
factor, defined as the signal detected in exact backscat- 
tering direction divided by the diffuse background, never 
exceeds two. This maximum value is reached if each pair 
of interfering paths has the same amplitude, and if single 
scattering can be suppressed. 

Concerning the nonlinear regime, previous studies have 
been restricted to the case of linear scatterers embedded 
in a uniform nonlinear medium 0, 0, 0- Here, it has 
been shown that the maximum enhancement factor re- 
mains two. As we will show in this letter, however, the 
situation drastically changes in the presence of nonlinear 
scattering (in contrast to nonlinear propagation) . In par- 
ticular, in the perturbative regime of at most one scatter- 
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FIG. 1: In the perturbative approach, we assume a single non- 
linear two-photon scattering event (□), but arbitrarily many 
linear scattering events (•). One of the two photons is fi- 
nally annihilated by the detector, thereby defining the pho- 
todetection signal, whereas the other one is scattered into an 
arbitrary direction. 



ing event with x^^^ nonlinearity, CBS arises from inter- 
ference between three amplitudes. Depending on the sign 
of the nonlinearity, this leads to an increase or decrease 
of the nonlinear CBS enhancement factor compared to 
the linear value two. The clearest manifestation of this 
effect, however, can be observed if the linear component 
is filtered out from the backscattering signal, such that 
the latter exclusively originates from nonlinear processes. 
Then, interference between three amplitudes allows val- 
ues of the CBS enhancement factor up to three. 

In the case of a disordered medium consisting of cold 
two-level atoms, the linear and nonlinear component of 
the backscattered light can be distinguished in terms of 
its frequency. This is due to the fact that scattering is 
purely elastic in the linear regime (described by single- 
photon scattering), whereas nonlinear multiphoton pro- 
cesses in general change the frequencies of the individual 
photons, giving rise to inelastic scattering |2i|. In this let- 
ter, we will restrict ourselves to the perturbative regime 
of weak laser intensities which is described by scattering 
of two photons 0- A typical scattering path is repre- 
sented in Fig. n Here, the two incoming photons prop- 
agate at first independently from each other to position 
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FIG. 2: In the presence of nonlinear scattering (□), there 
may be either (b) two, or (c) three interfering amplitudes con- 
tributing to enhanced backscattering, apart from single scat- 
tering (a), which only contributes to the background. In gen- 
eral, the case (c), which corresponds to maximum enhance- 
ment factor three, is realized if either both incoming photons, 
or one incoming and the outgoing detected photon exhibit at 
least one linear scattering event (•) besides the nonlinear one. 



r inside the disordered atomic medium, where they un- 
dergo a nonlinear scattering event. One of the two out- 
going photons then propagates back to the detector. The 
possibihty that the two photons meet again at another 
atom can be neglected in the case of a dilute medium, 
similar to recurrent scattering in the linear case 0. We 
can hence restrict our analysis to processes like the one 
shown in Fig. ^ with arbitrary number of linear scatter- 
ing events before and after the nonlinear one. As further 
discussed below, this perturbative approach is valid for 
small laser intensity and not too large optical thickness. 

To obtain the average intensity measured by the detec- 
tor placed in backscattering direction, we must identify 
those scattering paths whose interference survives the en- 
semble average over the random positions of the scatter- 
ers. Again, we consider the case of a dilute medium, 
where the typical distance between two scattering events 
(given by the linear mean free path £) is much larger 
than the wavelength A. First, we do not obtain inter- 
ference between diagrams where the nonlinear scattering 
event occurs at different atoms Just as in the lin- 
ear case, we do find interference between diagrams where 
the path of the detected photon is reversed. This gives 
rise to an enhancement of the average detection signal 
in exact backscattering direction, on top of the diffuse 
background. In contrast to the linear case, however, there 
may be up to three different interfering amplitudes, as ev- 
ident from the examples shown in Fig. [21 In general, the 
three-amplitudes case is realized if both incoming pho- 
tons, or one incoming and the outgoing detected photon 
exhibit at least one linear scattering event besides the 
nonlinear one. 

The perturbative calculation of the backscattering sig- 



nal relies on the fact that the scattering amplitude of a 
multiple scattering process like the one shown in Fig. ^ 
is obtained as a product of the scattering amplitudes for 
each individual scattering event. The single-photon scat- 
tering events can be treated by standard methods known 
from the theory of multiple scattering in the linear regime 
[lo|. which we shortly summarize as follows. Here, the 
central quantity is the average intensity /ij(r) inside the 
medium. It can be interpreted as the probability that an 
incoming photon with frequency w reaches position r af- 
ter arbitrarily many linear scattering events, and fulfills 
the following version of the radiative transfer equation: 

lUr) = e--/'+M\S^\^ [ dr'|G^(r,r')p/.(r'). (1) 
Jv 

Here, z denotes the distance from the boundary of the 
medium to r along the direction of the incident beam, and 
M the density of atoms which are randomly distributed 
in the volume V. Furthermore, 
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represents the scattering amplitude for a photon with 
frequency w = fee by a single atom, where F denotes the 
width of the atomic resonance at Wat- Eq. Q is valid in 
the near-resonant case of small detuning 6 ~ uj — Wat ^ 
Lu. For simplicity, we work with scalar photons, i.e., we 
neglect the vectorial nature of the light field. In this case, 
scattering is isotropic, with the total cross section given 
by a = 47r|iS'ijp. This is not a crucial assumption - the 
following treatment can be generalized to the vectorial 
case without serious difficulties. Finally, 
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describes average propagation between r and r' in the 
atomic medium, with refractive index given by 
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(4) 



The imaginary part of n^^ is determined by the linear 
mean free path £ — \/{Afa) at frequency uj. 

The first term in Eq. ^ represents the exponential at- 
tenuation of the coherent laser mode, i.e., light which has 
penetrated to position r without being scattered (Beer- 
Lambert law). The remaining term describes the diffuse 
intensity, i.e., light which has been scattered at least once 
before reaching r. The required solution of Eq. |^ is ob- 
tained numerically by iteration, starting from /^^(r) = 0. 
In terms of the average intensity the linear com- 

ponents of the backscattering signal are readily obtained 
by integration over the medium |lOl |. 

The second ingredient needed for our nonlinear pertur- 
bative analysis is scattering of two photons by a single 
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atom [Sj. Here, the scattered light contains an elastic 
component, with the same frequency w as the incident 
light, and an inelastic component with power spectrum 
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S' + iT/2 2S-d' + iT/2 
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where S' — to' — Wat denotes the final detuning. The ratio 
between the inelastic and elastic component is given by 
the saturation parameter s = 2r2^/(4<5^ + F^) (with the 
Rabi frequency Q), which measures the intensity of the 
incident laser in units of the saturation intensity 0|- 

Combining scattering of two photons by a single atom 
with the technique for treating the linear case outlined 
above, we arrive at the following expression for the inelas- 
tic background component at frequency uj' (postponing 
a more detailed derivation to a further publication): 

L,^{u;') = sPiu;') ^ (2/^(r) - er^^-/') I^,{v), (6) 

where A denotes the transverse area of the medium. The 
term 2/2(r) - e^^z/^ Q 

can be interpreted as the 
average squared intensity at position r inside the slab, 
which is larger than the square /^fr) of the average in- 
tensity due to speckle fluctuations The squared in- 
tensity induces a nonlinear atomic response, emitting a 
photon with frequency distribution Piuj'). Due to time 
reversal symmetry, the diffusion of this photon from r 
to the detector is given by the same expression I^^i (r) 
which describes diffusion of incoming photons to r. Note 
that, in the perturbative regime, the inelastic compo- 
nent contains nonlinear scattering [i.e. nonlinearity of 
the scattering cross section), but not nonlinear average 
propagation (i.e. nonlinearity of the mean free path). 
Since the latter leaves the photon frequency unchanged, 
it only affects the elastic component, which we assume 
to be filtered out from the detection signal. 

Concerning the calculation of the interference contri- 
bution giving rise to enhanced backscattering, the linear 
average intensity /cj(r), Eq. (Q, must be generalized in 
order to account for different frequencies uj and oj' in the 
interfering paths: 
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X / dr'G^(r,r')G:,(r,r')5c.,..'(r'). (7) 
Jv 

This function describes the ensemble-averaged product 
of two probability amplitudes, one representing an in- 
coming photon with frequency u propagating to position 
r, and the other one the complex conjugate of a pho- 
ton with frequency uj' propagating from r to the detec- 
tor. If w 7^ w', a nonvanishing phase difference between 
these amplitudes arises due to both scattering and aver- 
age propagation in the medium, since both the complex 
scattering amplitude, Eq. |(2Jl, and the refractive index. 



Eq. depend on frequency. In contrast, the phase 
difference due to free propagation [i.e. in the vacuum) 
can be neglected if F^ <C c, which is fulfilled for typi- 
cal experimental parameters [l^ . In the case w = w', 
5w,w(r) = Iuj[tc) reduces to the average intensity, see 
Eq. (Q. In terms of the iterative solution of Eq. iQ, we 
find the following expression for the interference term: 

a„(c.') = 4.P(c.') y^_[/^(r)|5^,^,(r)P 

-e-^/^Re{e»("--<')'=^5:,.'(r)} 



(8) 



with the linear mean free path at frequency lo' . In 
order to verify that the interference of either three or two 
amplitudes is correctly taken into account in Eq. ©, it 
is useful to write the average intensity / = exp(— C) + Ib 
as a sum of the coherent mode exp(— C) (with C = z/t) 
plus diffuse intensity Id- In the case w' = w, we obtain 
from Eqs. P5|) : 
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4/oe-^^ + 12/1)6-^+4/1)) , (9) 
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where the brackets denote the integral over the volume 
V of the medium, and (a,b,c) correspond to the three 
cases shown in Fig.|3 identified by different powers of dif- 
fuse or coherent light. As expected, the three-amplitudes 
case (c) implies an interference term twice as large as the 
background, whereas the case (b) leads to the factor 4/5 
known from the scalar two-atom solution (since one of the 
two interfering amplitudes is twice as large as the other 
one 0). Finally, as it should be, the single scattering 
term is absent in the interference term Cin. 

Using Eqs. (|6I8|I . we now calculate the CBS enhance- 
ment factor rj for the frequency component lo': 



ri[uj') = 1 + 



Cin(c^') 



(10) 
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The result is shown in Fig. 13 for a slab geometry with 
three different values of the optical thickness, b = 0.5, 1, 
and 2. Evidently, the largest values of the enhancement 
factor are obtained if the final frequency uj' approaches 
the initial one uj, since then the dephasing due to differ- 
ent frequencies vanishes. In this case, the enhancement 
factor is completely determined by the relative weights 
between the one-, two- and three- amplitudes cases shown 
in Fig. 12 cf. Eq. @ . As evident from the dashed line 
in Fig. 13 already at the rather moderate value b — 0.5 
of the optical thickness, the three-amplitudes case is suf- 
ficiently strong in order to increase the maximum en- 
hancement factor above the linear barrier rj = 2. With 
increasing optical thickness (and, if necessary, decreas- 
ing saturation parameter, in order to stay in the domain 
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FIG. 3: Spectral dependence of the enhancement factor for 
coherent backscattering from a slab with optical thickness h — 
0.5 (dashed line), & = 1 (solid line), and b — 2 (dotted line) in 
the perturbative regime of small saturation parameter s. The 
initial frequency ui = ojat is chosen on resonance (i.e. 5 = 0). 
The vertical line displays the frequency of elastically scattered 
photons, which must be filtered out in order to observe an 
enhancement factor larger than two. 

of validity of the perturbative approach, see below), the 
number of linear scattering events increases, which im- 
plies that the three-amplitudes case increasingly domi- 
nates, see Fig. |21 In this limit, the enhancement fac- 
tor approaches the maximum value three. At the same 
time, however, a larger number of scattering events also 
leads to stronger dephasing due to different frequencies, 
Lu' ^ uj. This results in a narrower shape of as a func- 
tion of uj' for larger optical thickness. Nevertheless, as ev- 
ident from Fig. 13 the enhancement factor remains larger 
than two in a significant range of frequencies lo' . Hence, 
an experimental demonstration of three-amplitudes in- 
terference requires a sufficiently narrow spectral filter, 
which should be placed as close to the laser frequency 
UJ as possible, but far enough to ffiter out the elastic 
component. We have checked that, without filtering, the 
enhancement factor decreases as a function of s (due to 
the negative sign of the elastic nonlinear component), in 
qualitative agreement with the experiment 12]. For a 
quantitative comparison, the elastic component, the ge- 
ometry of the medium, and the polarization of the pho- 
tons must be taken into account, which will be presented 
elsewhere. 

Finally, the domain of validity of our perturbative ap- 
proach, assuming a single nonlinear scattering event, re- 
mains to be discussed. A rough quantitative estimation 
can be given as follows: if pi (resp. P2+) denotes the 
probability for a backscattered photon to undergo one 
(resp. more than two) nonlinear scattering event, the 
perturbative condition reads p2+ <^ pi . If we assume for 
all scattering events the same probability s to be non- 
linear (thereby neglecting the inhomogeneity of the lo- 
cal intensity), we obtain pi ~ {N)s, and p2+ — (iV^)s^, 
where N denotes the total number of scattering events. 



and (...) the statistical average over all backscattering 
paths. Evidently, N and N'^ are expected to increase 
when increasing the optical thickness b. For a slab ge- 
ometry, we have found numerically that (N) oc b and 
(iV^) oc (in the limit of large b), concluding that the 
perturbative treatment is valid if sb^ <C 1. Let us note 
that we remain well below the threshold 2± 1 at which 
speckle fluctuations in a nonlinear medium become un- 
stable 13]. 

A natural way how to extend this work is to give up 
the perturbative assumption, and admit more than one 
nonlinear scattering event. This is necessary in order to 
describe media with large optical thickness, even at small 
saturation. Since the number of interfering amplitudes is 
expected to increase if more than two photons are con- 
nected by nonlinear scattering events, it may be possible 
to observe even larger enhancement factors - especially in 
the case of amplifying scatterers where the nonlinearity 
appears with positive sign. Finally, the effect demon- 
strated in this letter might also be relevant for strong 
localization in the nonlinear regime. 
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